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The determination of neutron star masses is reviewed in light of a new measurement of 
1.97 M© for PSR J1614-2230 and an estimate of 2.4 M© for the black widow pulsar. Using 
a simple analytic model related to the so-called maximally compact equation of state, 
model-independent upper limits to thermodynamic properties in neutron stars, such as 
energy density, pressure, baryon number density and chemical potential, are established 
which depend upon the neutron star maximum mass. Using the largest well-measured 
neutron star mass, 1.97 Mq, it is possible to show that the energy density can never 
exceed about 2 GeV, the pressure about 1.3 GeV, and the baryon chemical potential 
about 2.1 GeV. Further, if quark matter comprises a significant component of neutron 
star cores, these limits are reduced to 1.3 GeV, 0.9 GeV, and 1.5 GeV, respectively. 
We also find the maximum binding energy of any neutron star is about 25% of the rest 
mass. Neutron matter properties and astrophysical constraints additionally imply an 
upper limit to the neutron star maximum mass of about 2.4 Mq. A measured mass of 
2.4 Mq would be incompatible with hybrid star models containing significant proportions 
of exotica in the form of hyperons, Bose condensates or quark matter. 

Keywords: neutron star masses; equation of state. 

1. Introduction 

One of the most fascinating stories in astrophysics concerns the accumulation of 
precisely measured neutron star masses. Gerry has had a long-time interest in these 
measurements, as for many years he has maintained that the neutron star maximum 
mass is no greater than (1.5-1.6) Mq. From a theoretical perspective, a maximum 
mass of 1.5 fits his fondness for the effects of kaon condensation proposed by 
Kaplan and Nelson HI the effective kaon mass falls with increasing density and the 
eventual onset of kaon condensation at a few times the nuclear saturation density 
softens the equation of state (hereafter, EOS) and leads to a rather small maximum 
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mass. He takes as observational support of this thesis the "missing" neutron star in 
the remnant of SN 1987A, which had a baryon mass not in excess of 1.7 Mq based 
on estimates of the initial mass of the star and its ejected mass He reasons that 
the maximum stable mass of the proto-neutron star as it deleptonized must have 
eventually fallen below its actual mass, prompting a collapse to a black holJ^l. Of 
course, this had to have happened more than 12 seconds after the core collapse since 
neutrinos were observed during that time period and their emission is believed to 
abruptly cease 3^ once the black hole's event horizon forms. 



Measurements of Neutron Star Masses 

The most accurate measurements of neutron star masses are for pulsars in bound 
binary systems. In these systems, five Keplerian parameters can be precisely mea- 
sured by pulse-timing techniques including the binary period P, the projection 
of the pulsar's semimajor axis on the line of sight Op sinz (where i is the binary in- 
clination angle), the eccentricity e, and the time and longitude of periastron Tq and 
uj. Combining two of the observational parameters, one can form the mass function: 

/27ry/apsiniy Mo (M^nif 

where M = Mp + Mc is the total mass, Mp is the pulsar mass, and Mc is the 
companion mass (all measured in Mq units). The constant Tq = GMq/c^ is 4.9255 
11 s. The mass function fp is also equal to the minimum possible mass Mc for the 
companion. 

The inclination angle i is often the most difficult aspect to infer, but even if it 
was known a priori the above equation only specifies a relation between Mp and 
Mc unless the mass function fc of the companion is also measurable. This occurs 
in the rare case when the companion is itself a pulsar or a star with an observable 
spectrum. 

Fortunately, binary pulsars are compact systems and general relativistic effects 
can often be observed. These include the advance of the periastron of the orbit 

^ = ^{^y\MTQ)^/Hl~eT\ (2) 

the combined effect of variations in the tranverse Doppler shift and gravitational 
redshift around an elliptical orbit 

f Py^' Mc{M + Mc) ^2/3 

^ = ^2^; MV3 , (3) 

and the orbital period decay due to the emission of gravitational radiation 

P (^_lE^y'' (, + 73^. ^ 37 \ _ , y.MpMc 

5 [ P 96 r ' Ml/3- 
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The inclination angle can be constrained by measurements of two or more of these 
effects. However, only in extremely compact systems is this precisely possible. Oth- 
erwise, additional effects, such as an eclipse or limits obtained from the lack of 
an eclipse, or Shapiro time delay, which is caused by the propagation of the pul- 
sar signal through the gravitational field of its companion, must be obser ved. The 
Shapiro ^ delay in general relativity produces a delay in pulse arrival time 



Ssicj))^2M,TQ In 



1 + e cos ( 



1 — sin(w + ( 



(5) 



where (f) is the true anomaly, the angular parameter defining the position of the 
pulsar in its orbit relative to the periastron. 5s is a periodic function of with an 
approximate amplitude, 



In 



1 + e sin uj\ / 1 + sin i 



1 — e sin w / V 1 — sin i 



(6) 



which is large only for edge-on binaries, sini ~ 1 or those which have both large 
eccentricities and large magnitudes of sinw. Only a fraction of pulsars in binaries 
have two or more sufficiently well-measured relativistic effects to enable precise 
measurements of the pulsar mass Mp . 

In some cases, the companion star can be optically detected. In most cases, the 
companion is a white dwarf. It is possible to estimate white dwarf masses from 
observations of the optical flux, effective temperature and distance, with the latter 
being especially problematic. Another possibility is to measure the surface gravity 
from spectral measurements, but currently this is subject to large systematic effects. 

Less accurate measurements of neutron star masses are possible in X-ray bina- 
ries, in which X-rays are emitted by matter accreting onto a neutron star from a 
companion star. Both X-ray and optical observations can yield both mass functions 
fp and /c, but fp is subject to large uncertainties due to the faintness of optical ra- 
diation. In the case of the black widow pulsar, PSR B1957-I-20, optical observations 
of the companion yield both a mass function and an estimate of the inclination i 
from the shape of the light curve The inferred semimajor axis of the companion's 
orbit has to be corrected for the finite size of the companion, however (see below), 
introducing a source of systematic error in the pulsar mass determination o f th is 
system. Nevertheless, interesting lower limits to the pulsar mass are obtained 



Recent Measurements 

The most recent information, as of November 2010, is summarized in Figure [Hand 
Table [Hand collects information from Refs.fTOltolSOl This compilation represents a 
significant update to the figure of observed masses and references in Ref. [52] 

There is now ample observational support from pulsars for neutron stars with 
masses significantly greater than 1.5 Mq. These include PSR J1903-f 0327 which has 
a main-sequence companion; the pulsars I and J in the globular cluster Ter 5, the 
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Fig. 1. Measured neutron star masses. References in parenthesis following source numbers are 
identified in Table [Tl 
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Table 1. Neutron star mass measurements with Icr uncertainties. Reference letters correspond to Fig. 
1. For each group, mean and weighted masses are indicated. 
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globular cluster binaries PSR J1748-2021B and PSR B1516-I-02B, and PSR J1614- 
2230, all with white-dwarf companions; and the black widow pulsar B1957-I-20. The 
inclination angles of the two binaries containing pulsars I and J in the globular clus- 
ter Ter 5 are unconstrained by observation, but if their inclinations are individually 
randomly directed, there is a 95% chance that at least one of these pulsars is greater 
than 1.68 Mq ^^SI. Assumptions related to random orientation also have roles in the 
mass determinations for PSR J1748-2021B and PSR B1516-I-02B. One has to be 
careful about the assumption of random orientation, however, as the selection of 
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potential candidates for follow-up observations necessary for accurate mass deter- 
minations are biased. The case of PSR J0751+1807 illustrates the point: initially 
estimated to have a mass 2.2 ± 0.2 M©, the mass was recently established to be 
1.26 Mq upon measurement of another relativistic parameter. Nevertheless, about 
2 years ago, a measurement of the inclination of PSR J1903-f 0327132! using Shapiro 
time delay lead to an accurately determined mass of 1.67±0.01 M0, replacing PSR 
B1913-I-16 as establishing the value for the minimum maximum mass. 

But even more outstanding is the recent determination of 1.97±0.04 Mq for the 
mass of PSR J1614-223c|3SI, also using Shapiro time delay to measure the inclination. 
This 3.15 ms pulsar is in an 8.7 d nearly circular orbit about an 0.5 Mq companion, 
with CpSini = 11.3 light-second and sini = 0.99989, i.e., it is almost edge on. The 
Shapiro time delay amplitude, from Eq. (|6]), becomes 

As ~ 4:McTq ln(2/ cos i) ~ 48.8 /xs . (7) 

By virtue of its accuracy, this mass measurement has become the new standard for 
the minimum value of the neutron star maximum mass. 

In addition, a number of X-ray binaries seem to contain high-mass neutron 
stars: about 1.9 Mq in the case of Vela X-1 and 2.4 Mq in the case of 4U 1700-377. 
Nonetheless, the large systematic errors inherent in X-ray binary mass measure- 
ments warrant caution. 

The black widow pulsar represents an intriguing case, with a best estimate of 
about 2.4 Mq. This system has both pulsar timing and optical light curve infor- 
mation. It consists of a 1.6 ms pulsar in a nearly circular 9.17 h orbit around an 
extremely low mass companion, Mc — 0.03 Mq. The pulsar is eclipsed for about 
50-60 minutes of each orbit, but considering that OpSini = 0.089 light-second = 
0.038 Rq, and OcSini ~ 3 Rq is Mp/Mc — 80 times larger, the eclipsing object has 
to be much larger (~ 3 Rq/10) than the size of the companion star. It is believed 
that irradiation of the companion by the pulsar strongly heats its nearside to the 
point of ablation, leading to a comet-like tail and a large cloud of plasma which 
is responsible for the eclipsing. The pulsar is literally consuming its companion, 
hence the name "black widow" , and has reduced its mass to a small fraction of its 
original mass. The irradiation also produces an enormous (factor of 100) variation 
in the brightness of the companion depending on how much of its nearside is visi- 
ble during its orbit. The companion is bloated and nearly fills its Roche lobe. The 
companion's optical light curve allows the mass ratio Mp/Mc and the inclination 
angle i to be estimated. However, the large size of the companion means that the 
"center of light" of the system is not equivalent to its "center of mass" : the optical 
light curve depends on the projected semi-major axis of the irradiated nearside of 
the companion rather than the projected semi-major axis of the center of mass of 
the companion. The extreme cases are either that the companion has zero radius or 
that it completely fills its Roche lobe, but estimates based on modeling consider- 
ably reduce the allowed range. The extremes lead to a range 1.7 < Mp/MQ < 3.2, 
but the likely value is 2.4 ± 0.4 Mq. It will clearly be valuable to extend observa- 
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tions and modeling of this system since a 2.4 Mq neutron star would have profound 
implications. 



What Gerry Would Say 

We would love to be able to argue with Gerry about these observations and their 
implications, and so we felt this was an ideal topic to explore in this contribution. 
We can, however, guess what Gerry would have said: "We can't tell nature how 
to behave; we have to explain it" and, as he often quoted Bismarck, "Was schert 
mich mein Geschwaetz von gestern?" ("Why do I care about my twaddle from 
yesterday?" ) 

What We Say 

One of the most surprising consequences of a large pulsar mass measurement is 
that it severely restricts values of the maximum central density and pressure in 
that star, and, by extension, that of any other neutron star as well. This realization 
was elaborated in Ref. [52l but it can be extended to restrict the possibilities for the 
presence of exotic matter, such as hyperons, Bose condensates and/or quark matter, 
in neutron star interiors. Moreover, we show here that upper limits to the chemical 
potentials of quark or baryonic matter can be established using the available mass 
data, and that these are smaller than often employed in calculations of the properties 
of quark matter. The consequences of well-determined large neutron star masses is 
the topic to which we now turn our attention. 



2. Neutron Star Structure 

The global properties of neutron stars are determined, ass uming general relativity, 
by the Tolman-Oppenheimer-Volkov (TOV) equations 1^ 

dp _ G{e + p){m + Anr^p/c^) 

d?^ r(rc2 - 2Gm) ' ^' 

^ = We/c^ , (9) 



dN , / 2Gm\ , , 

= 1 , (10) 



-1/2 

dr \ r(? ) 

where m is the enclosed gravitational mass, M is the enclosed baryon number, p is 
the pressure, e is the total energy density, n is the baryon number density, and r is 
the radial variable. Choosing a value of the pressure or baryon or energy density at 
the star's center, and integrating to the surface where the pressure vanishes, results 
in a specific mass and radius, as shown in Fig.[2j A consequence of general relativity 
is the existence of a maximum mass: configurations with larger central pressures are 
dynamically unstable and collapse into black holes. 
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Fig. 2. Solution of the TOV equations. Integrations beginning at the star's center with pressures 
and baryon number densities (in units of the nuclear saturation density Us) shown by the dots 
in the left panel result in specific masses and radii, shown by the corresponding dots in the right 
panel. 



It is convenient to recast the TOV equations in terms of the dimensionless vari- 
able h defined hj dh = dp/{e + p) = dn/fj. where fi = {e + p)/n is the baryon 
chemical potential, which in beta equilibrium is the same as the neutron chemical 
potential. Therefore, 

fj, = Hoe^ , (11) 

where yUo = So/no, and Hq are the chemical potential, energy density and baryon 
number density, respectively, where p = h = 0. For a star with a normal crust, 
fio — 930 MeV, the approximate mass-energy of iron nuclei per baryon. Note that 
the quantity h is related to the metric function v, where gtt = and dh = —du/2. 
One finds 

dr'^ _ 2r'^c^ 
'dh ~ G~ 
dm A-KT^e 
'dh ^ G~ 

dAf _ 47rr^/^ 
dh ~ G 

Writing the TOV equations this way ensures that near both the center and the 
surface the right-hand sides are finite. The equations are integrated from a central 
value h = he to h = 0. In order to obtain a specific total gravitational mass M = 
m{h = 0), one iterates to find the appropriate value of he- For a given EOS, that is, 
the functions e{h) and p{h), a specific value or he, hmax, results in the maximum 
mass configuration where M = Mmax, N = N„iax and R = R„iax where R = r{h = 
0) is the stellar radius and A^" = M{h = 0) is the total baryon number. The baryon 



- 2Gm 
+ Anr^p 
■ - 2Gm 



nc 



\/r<P' — 2Gm 
mc^ + Anr^p 



(12) 
(13) 
(14) 
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Fig. 3. The solid line is a typical neutron star EOS. The dashed line shows the maximally compact 
EOS formed by minimizing the pressure p at low energy density e (i.e., p = 0) and maximizing 
p, subject to causality, at high density (i.e., p = e — So), p and e are in the same arbitrary units, 
with So in this example set to 2. 

number is useful in the definition of the binding energy 

BE = {Nrrib - M)c^ , (15) 

where is the baryon mass. 

It is an inescapable consequence of general relativity that all EOSs result in a 
maximum mass configuration, in contrast to Newtonian gravity for which no such 
limit exists. It is intuitive that the largest maximum masses will be a consequence 
of the stiff est EOSs, i.e., those with the largest pressures for a given density. It is 
commonly assumed that the maximum realistic limit to the pressure is obtained by 
assuming causality, in other words, that the sound speed is limited everywhere by 
light speed. The assumption of causality led Rhoades and Ruffini to their estimate 
for the maximum neutron star mass. It is logical, therefore, that the maximum mass 
configuration so obtained must also have the largest possible values for the central 
pressure, energy density and chemical potential. The main motivation of this paper 
is to explore the consequences of this in light of new observations of large neutron 
star masses. 

It was argued by Haensel and Zdunik^SSland by Koranda, Stergioulas and Fried- 
man ESI that the most compact configurations result from the combination of a soft 
EOS at small densities and a stiff EOS at large densities. The extreme limit of this 
situation is shown in Fig. |3l where the pressure vanishes below the density £o and 
is assumed causal at all higher densities (the so-called maximally compact EOS): 

p = £ < £oi p — e ~ £o £> £o- (16) 

Since £o is the only EOS parameter, the TOV equations scale with it according 
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to 153 

-1/2 -1/2 /ir,\ 

e oc Eo, p (X Eo, n (X So, m (x ' , r oc ' . (17) 

Following Lattimer, Prakash, Masak and YahilEU (hereafter LPMY), we wish to 
study a somewhat more general EOS, which encompasses the maximally compact 
EOS as a special case, namely, 

p = s{e-eo), (18) 

where s is a constant representing the square of the adiabatic sound speed. We can 
then express e, p and the baryon number density n in terms of h as 

£ = t-ie„(s-ie''* + 1) , p^t'^Eoie''* -1), n ^ Eoc''^' / fio (19) 

where t = I + 1/s. Using the dimensionless variables 

x = r^eoG/c\ y = mt^J^G^^ I c\ z = Nel'^ HoG^'^ /c"^ (20) 

the scale-free TOV equations become 

dx 
dh 
dy_ 
dh 
dz 
dh 

For a given value of s, the value of hmax is the value of he that maximizes the 
total gravitational mass y{h = 0). By virtue of Eq. (ITU) . h„iax then determines 
the maximum possible chemical potential for that EOS. Table [5] summarizes the 
properties of the maximum mass solutions for the cases s = 1, 4/5, 2/3 and s — 1/3. 
The case s — 1 represents the maximally compact star; as discussed in §4, the case 
s = 1/3 is potentially applicable to quark stars or stars with quark cores. For each 
s are listed the central values of h, e and p, and the surface values of a:, y and z, 
together with other dimensionless quantities of interest. 

We will demonstrate that the maximally compact and the s — 1/3 EOSs result 
in useful approximate limits for the maximum values of thermodynamic quantities 
in the interior of neutron and quark stars as well as their maximum masses. 

3. Maximally Compact EOS 

We first examine the maximally compact EOS, i.e., the case s = 1. Expressing 
the results for s = 1 shown in Table [2] in astrophysical units (note that 1 MeV 



2x(xi/2 - 2y) 
y + Airx^/'^p/eo 

So dh 



= 27r 



dx 



So (x^/^ -2yy/'^ dh' 



(21) 
(22) 
(23) 
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= 1.323 X 10"*^ km-^ and 1 Mq 1.477 km), we find 
Mma. = yma^c^e^'/'G-^/s ^ 4.09 (e./eo)!/^ Mq , 

Rma. = Xraa^c^t-^I^G-^I^ = 17.07 (e./so)'/' km, 

£max = (e^''-- + l)£o/2 - 3.034 £o , 
Pmax = (e^''-^^^ - 1)£„ = 2.034 , 
f-^max — /^o^ — 2.251 /io , 

?lmax = {Smax + Praax)e~^""°-" / i^o = 2.251 (Eo/Mo) • (24) 

where Eg — TO&f^s — 150 MeV is the mass-energy density at normal nuclear satu- 
ration baryon density ~ 0.16 fm^'^. The result for Mmax is very close to, but 
slightly smaller than, that discussed by Rhoades & Ruffini the difference being 
their substitution of a low-density neutron star crust EOS for densities less than 
£/, where Sf > £o and pf > 0. However, the difference is very small, showing the 
validity of the maximally compact EOS for an approximate limit. 



s 


hfnax 


^max / 


Pmax / 




Umax 






^raax 1 Umax 


1 


0.8116 


3.034 


2.034 


0.05779 


0.08513 


0.1127 


0.02199 


2.824 




0.855 


3.00 


2.00 


0.0635 


0.0854 


0.109 


0.0219 


2.95 


4/5 


0.727 


3.29 


1.83 


0.0535 


0.0787 


0.103 


0.0204 


2.94 


2/3 


0.663 


3.55 


1.70 


0.0500 


0.0732 


0.0946 


0.0190 


3.05 


1/2 


0.569 


4.01 


1.51 


0.0443 


0.0644 


0.0810 


0.0166 


3.27 


1/3 


0.4520 


4.826 


1.275 


0.03649 


0.05169 


0.06212 


0.01289 


3.696 


1/3* 


0.4350 


6.585 


1.639 


0.02696 


0.04307 


0.05134 


0.01221 


3.813 


1/3** 


0.491 


4.68 


1.23 


0.0439 


0.0522 


0.0607 


0.0128 


4.01 



The results for Emax/so 8.nd pmax/so for the case s = 1 were previously obtained 
by Ref. [56l What they did not apparently realize was the significance of the quan- 
tity EmaxUmax/^o, which, when expressed in astrophysical units for the maximally 
compact case s = 1, is 

£,naxM^^^ = 1.358 X 10^6 g ^2^ (25) 

This equation, of course, states that the central energy density of the maximum 
mass star is inversely proportional to the square of the maximum mass. Since the 
maximally compact EOS very closely predicts the largest value for the maximum 
mass, it seems to follow that this equation should also describe the maximum pos- 
sible energy density as well. In other words, substituting the mass of the largest 
observed neutron star mass for Mmax in the above expression not only determines 
the largest energy density in that star, but it might yield an approximate upper 
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Approximate Central Baryon Density n /n 
012345 6 7 8 9 10 11 12 13 14 




1 2 3 4 5 6 

Centroi Moss-Energy Density (10^^ g cm"^) 

Fig. 4. The central density versus mass relations for the maximally compact EOS {s = 1, Eq. I|25|l ') 
in comparison to Tolman VII + causality and the s = 1/3 (Eq. I I39I1 ) self-bound cases. Various 
points are results from different kinds of EOSs, as labelled. Horizontal dashed lines indicate masses 
of 1.4, 2.0 and 2.4 M©, respectively. 

limit to the energy density to be found in any neutron star. Should a more massive 
star ever be measured, the upper energy density limit would have to decrease. 

It is interesting that the analytic solution to the TOV equations for the energy 
density distribution 

e^e,[l-{r/R)% (26) 

i.e., the Tolman VII solutionis, coupled with the causality condition Eq. (|27|) . 
predicts a value for £maxM^^^ to within a few percent of Eq. (pS)) . even though 
the EOS resulting from Eq. ([26|) becomes acausal at sufficiently large densities. 
Ref. [SH showed that the limit established from the Tolman VII solution is satisfied 
phenomenologically by all realistic EOSs, as illustrated in Fig. HI By comparison, 
the limit defined by Eq. (|25|) is slightly violated by two of the EOS models but 
nevertheless represents a realistic limit. As was the case for the estimate of the 
maximum neutron star mass, the slight violations result from the assumption of 
zero pressure at low densities. 

Ref. [56] argues that the maximally compact EOS generally predicts the smallest 
radius for a configuration with arbitrary mass. In this event, the minimum com- 
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pactness R/AI for any neutron star is determined by 

RmajMraax = 2.82AG/c' (27) 

a result consistent with discussions by Refs. [5^ [551 [551 and IBOl Modifying the EOS 
at low densities, i.e., for e < (2 — 3) Ss, does not significantly alter the values for 
Mmax, but the values for Rmax are moderately increased!^. 

Ref. [Sni also demonstrated that the axially symmetric structure equations in 
general relativity also scale exactly for the maximally compact EOS in terms of 
Co- The maximum rotation rate, limited by mass-shedding at the equator, scales as 
^y Mmax I Rfnax ^ where the mass and radius refer to the values for maximum 
mass of the non-rotating configuration. As a result, they calculated the minimum 
spin period for a causally-limited EOS: 

\ ^ f Rmax \ ^ Mmax 



P™„ = 0.74 -p^ nis - 0.20^:^ ms (28) 

Lattimer and Prakash^^ showed that for realistic equations of state the coefficient 
0.74 ms is effectively replaced by 1.00 ms, a value almost exactly equal to the 
general relativistic result in the case that the gravitational potential of the rotating 



star is due to a point mass at the origin Furthermore, the revised formula 
accurately predicts the maximum rotation rate of any star, not just the maximum 
mass configuration, if M^ax and Rmax are replaced by their non-rotating values. 

The maximally compact configuration is also expected to maximize the binding 
energy, 

^max Vrnax 1 TTo ~ — 0.337 ]\^IfyiaxC • V^^) 
Mo / £y'G3/2 

The maximum binding energy is thus 25.2% of the rest mass, a fraction that is 
independent of Eq and Mmax- 

It is intersting that each of the cases shown in Table [2] yields limiting values for 
hmax and therefore fj-max that are also independent of Eq and therefore also M^ax ■ 
In the case of the maximally compact EOS, one finds the upper limit to the baryon 
chemical potential 

fimax = 2.251 X 930 MeV = 2093 MeV . (30) 

These limits to the chemical potential, baryon number density, and binding 
energy have not been previously calculated. We next explore the implications of the 
large neutron star masses recently measured, and whether the limits obtained from 
the ideal s = 1 and s = 1/3 cases pose realistic upper bounds for thermodynamic 
quantities within neutron stars. 

Implications of the newly-measured mass 

The measured mass of 1.97 Mq can now be used to set specific limits because the 
value of eo is limited according to the first of Eqs. ([M)) : < 4.33 Eg- This implies 
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that 

emax < 3.034 X 4.33 = 13.1 Ss = 1.97 GcV fm^^ , 
Pmax < 2.034 X 4.33 Es = 8.81 Es = 1.32 GcV fm"^ , 
n,nax < 2.251 X 4.33 Us = 9.75 1.56 fm"^ . (31) 

Recall that the upper hmit to the chemical potential, Umax — 2.23es/"-s = 2.1 GeV, 
is independent of the measured mass. 

Finally, let us note the implications of a hypothetical measurement of 2.4 M0, 
for example, were the current estimate of the black widow pulsar to be confirmed. 
In this case, by virtue of the first of Eqs. we find So < 2.92 Ss- This in turn 

implies that 

emax < 3.034 X 2.92 = 8.86 = 1.33 GcV fm"^ , 
Pmax < 2.034 X 2.92 = 5.94 = 0.89 GcV fm"^ , 
Umax < 2.251 X 2.92 Us = 6.57 n, = 1.05 fm"^ . (32) 

4. Quark Matter EOSs 

As at sufficiently high density a quark-hadron transition is likely to occur, it is 
reasonable to consider a quark matter EOS for the highest densities that might 
exist in neutron stars. In the case that quark matter has a lower energy per baryon 
at zero pressure than iron nuclei, a pure quark star is possible in which case the 
pressure vanishes at non-zero baryon density . This is an example of a self-bound 
star. The radii of these stars decrease with decreasing mass in contrast to the case 
of normal neutron stars. The existence of a thin hadronic crust on such a star is 
possible, but this would not affect considerations of the maximum mass or central 
density, pressure or chemical potential of the star, although the mass-radius relation 
for small masses is dramatically altered. 

The term hybrid star refers to the case in which the star's interior contains 
exotica such as hyperons, Bose condensates and /or quark matter whereas the star's 
exterior is made up of hadronic matter as in a normal neutron star. The presence 
of matter extending to vanishing density causes the radii of hybrid stars to increase 
with decreasing mass as for normal neutron stars. A quark matter EOS is obviously 
softer than the maximally compact EOS, and since the maximum mass depends 
primarily on the EOS at the highest densities, the maximal constraints for pure 
quark stars or for neutron stars with quark cores will be more severe than for the 
maximally compact case. A relevant question is whether or not the transition to 
quark matter is at low enough density to occur in a neutron star. 

MIT Bag Model 

Quark matter is likely to be highly interacting for densities near the transition 
density. Nevertheless, the leading order interaction terms are generally attractive 
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SO that a reasonable upper limit to the pressure of quark matter is that for non- 
interacting massless quarks, 

p=i(e-4B), (33) 

where B is the MIT bag constant (see, e.g., Ref.|Ml)- This is just the EOS of Eq. ^ 
with s = 1/3 and Co — 4i?. Even if QCD corrections are incorporated, the generic 
form in Eq. (IT51) is not greatly altered Additional corrections for a non-zero 
strange quark mass and for pairing introduces attractive contributions and must 
lower the maximum mass so we ignore these for now. If the EOS of quark matter 
approaches the ideal of Eq. (1551) at high density, and there is significant quark matter 
in a neutron star interior, we conjecture that limits observed in the case s = 1/3 
could more realistically be expected to approximately constrain the maximum mass 
and maximum possible energy density, pressure and chemical potential than what 
is determined from the maximally compact EOS. Irrespective of the value for s, 
the relation h — \n{^/^o) is still valid. The limits to mass and compactness were 
investigated for the case s = 1/3 by LPMyI^. The results for the dimensionless 
integrations are contained in Table [2j All results for structural variables scale with 
So just as in the maximally compact case. 

^From Table [2] it is clear for the case s = 1/3 that the maximum value of h, 
hmax = 0.4520, is considerably smaller than its value for the maximally compact 
EOS. The maximum mass stars are also not as compact, the limiting R/M being 
3.696. The various dimensionless limits for the case s = 1/3 include 

M„,, = 2.48 {es/£of'^ Mq, 
N^axTTib = 2.98 {mbC^/^jio){es/eoY''^ , 
Rmax = 13.56 [Es/eof^ km, 

fJ-max = /ioc''"""" — 1.571 /io 7 

Emax = (e^''-- + l)eo/2 = 4.826 , 
Pmax = (e^''""- - l)eo = 1.275 , 

^max — i^max ~^ Pmax^ / f-^max — 3.884 {^Sq / l^o) • (^4) 

The results for Mmax, Rmax and Smax were previously obtained by Witten^SZl^ ^ 
result of great significance is that, under the reasonable assumption that Eq > Eg 
(otherwise deconfined quark matter would be visible in atomic nuclei) , the maximum 
mass of a pure quark star or a hybrid star with a substantial quark core can be no 
greater than 2.48 Mq. We will see that, in the case of pure quark stars, more realistic 
quark models generally have smaller pressures for a given energy density than the 
MIT bag model and therefore smaller maximum masses. In the case of hybrid stars, 
if quark matter exists in the cores, the limiting EOS is essentially that of pure quark 
matter. The magnitude of the transition density rit depends sensitively on whether 
the hadronic EOS is soft or stiff at the relevant densities. For stiff hadronic EOSs fit 
is often larger than the central density of the maximum mass star and quarks will 
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not appear. In the case that nj is less than several times n^, we will show that the 
hadronic EOS at lower densities does not substantially alter the maximum mass. 
Thus, the observation of a 2.4 Mq neutron star, with the certainty that the true 
maximum mass is larger, will greatly challenge the notion of the existence of quark 
matter in neutron star interiors. 

It is important to note that the chemical potential fi = /Xq In refers to the 
baryon chemical potential and n refers to the number density of baryons. For neutral 
quark matter with three flavors, the quark chemical potential is 1/3 the baryon 
(neutron) chemical potential and the quark density is 3 times the baryon density. 
Nevertheless, the value of nfi is the same for both types of matter. The value of Ho 
will depend on assumptions regarding the surface of quark matter stars. It is widely 
believed that a thin baryonic crust will exist, and in this case, the value of fj,o will 
be the same as previously considered: Iron at zero pressure with /io ~ 930 MeV. If 
Ho = 930 MeV, then we immediately find for the case s = 1/3 that fimax = 1-46 
GeV. As we show below, this limit for a self-bound MIT bag quark star remains 
approximately valid for self-bound interacting quark stars and for hybrid quark 
stars with a hadronic exterior. 

Importantly, the TOV equations do not depend on the value for /io as can be seen 
from Eq. ()13p . but absolute values for the number density and chemical potential 
do depend on this quantity. In order to discuss separate limits for the chemical 
potential or number density, one has to have a specific model for the dependence of 
the energy density on the chemical potential. For massless quarks in the MIT bag 
model with bag constant B, the relation is 

where — ^/i is the quark chemical potential. Each of the three flavors of quarks 
has the same chemical potential and density in the massless, neutral, case. We can 
therefore directly use the results from Table [5] using — 4-B to obtain 




Note that the limiting value for /igo < 310 MeV fom the condition of stability of 
quark matter at zero pressure implies B^^'^ < 163 MeV for the MIT bag model 
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Constraints from the 2-solar-mass star 

We now see that the observed mass of 1.97 M© suggests that 

So < (2.48/1.97)2e, = 1.59 e^, B^/^ < 146.3 MeV 

No< (^^^^^^y\hcf/' = 278.7 MeV, 

1/4 /I r;o=- \3/4 

= 0.856 fm"^ (37) 



fl.59e 

,r2 



he 

Therefore, the MIT bag model (without strong interaction corrections) cannot ex- 
plain the 1.97 Mq star unless strange quark matter is stable. The condition that 
fiqo > 310 MeV, necessary for iron to be the ground state of matter at zero pressure, 
implies that the maximum mass is 1.58 M© if the MIT bag model is correct. 
One also finds 

Sma^ = 1.59 X 4.826 e, = 7.669 e, , 
Pmax = 1.59 X 1.275 £, = 2.026 e, , 

t^q,max = 1.571/iqo = 437.8 McV , 

Pmax H" ^max ^ r\nn /oo\ 

nb.niax ^ — ^ = 6.926 Us , (38) 

and 

SmaxMl^,^ = 7.91 X 10^5 g cm-3 M| . (39) 

This relation is shown in Fig.|4l The 1.97 M© star thus limits the energy density to 
2. Ox IQi^ g cm~3 = 1.13 GeV fm~'^, just about half the value for the maximally com- 
pact EOS. Likewise, the maximum neutron chemical potential is ^max — '^fJ-q.max 
or 1.29 GeV which is considerably less than the absolute upper limit of 2.09 GeV 
from causality alone. 

Note that the reasonable constraint Eo > Ss implies B^^^ > 130.3 MeV in this 
model. The relatively narrow restrictions on the value of B imply that 248.2 MeV < 
/igo < 278.7 MeV. If one further assumes that the baryon density at the surface 
of the star is at q times the nuclear saturation density, or = Sqris, one has 
1.26 < q < 1.78. If a neutron star of 2.4 M© were to be found, this range would 
shrink to an extreme degree. 

Modifications to the Bag Model 

In the model developed by Alford et al. ESI, QCD corrections are included with a 
parameter , and the effects of a finite strange quark mass and pairing are included 
with a term proportional to the parameter a2- 

where Beff is an effective bag constant. 
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We first consider the case when 02 =0. Obviously, the introduction of the pa- 
rameter 04 and the replacement of B with B^f / will not affect the structure of the 
star and the dimensionless results from the case s = 1/3 still apply as for the MIT 
bag model. The introduction of the parameter 04 affects the chemical potential and 
number densities, however. The quark chemical potential and quark number density 
at the surface of the star where p = can now be determined for 02 = (see also 
Ref.ESl): 



3 



(±iy^^^ 0.766 --^^ aT fm-^(41) 
\hc) \ 141 MeV / ^ ^ ' 



In this case the observed mass of 1.97 M0 suggests that 

Hqo < [^"^^y = 278.7 a-'/' MeV, 



304 ^L59e,V/^ _ , m 



he 



3 



"9° < 1 7^ ) ( -7-^ ) = 0-856 04 fm" ■ (42) 



It is now evident that 04 < 0.656 is required if both fiqa — 310 MeV and M^ax > 
1.97 M0 in this model. Further, one has 

tiq.,max = 1.57lAiqo = 437.8 al^'^ MeV, 

Pmax ~^ £max „ 1/4 / a o\ 

nh,max = — ^ = 6.926 a/ Us . (43) 

The introduction of the 04 term does not alter the hmiting energy density inferred 
from the observed 1.97 Mq star. However, the maximum baryon chemical potential 
1.29/ay^ GeV increases moderately from the noninteracting quark case. 
Nevertheless, the condition that jiqo = 310 MeV does not change the limiting value 
of Umax — 1.46 GeV from that case. We note that one expects from theories of 
interacting quark matter that 04 0.7 

The constraint £0 > Ss implies -B^^^ > 130.3 MeV in this model, unchanged by 
the a4 parameter. The restriction on the value of B^ff from the observed 1.97 Mq 
star now implies that 248.2 < fiqoO,^^^ < 278.7 MeV. Finally, assuming the baryon 
density at the surface of the star to be qus, one finds 1.26 < qa_^ ^^"^ < 1.78. 

We next have to determine how modifying the quark EOS by introducing the 
02 term changes these results. The 02 term is generally a small correction to the 
quark EOS. In solving the TOV equation, we may once again scale energy and 
pressure with Bo — 'iBeff, but there is an additional parameter 02. Furthermore, 
So is no longer the value of e when p ^ 0. So the introduction of 02 is a moderate 
complication. 
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Although there are now 3 parameters affecting the quark EOS, there is a degen- 
eracy in the parameters as far as structural considerations are concerned. We write 
Eq. (1401) in the dimensionless form 



B. 



B. 



eff 




(44) 



Therefore it is only necessary to consider as parameters the quantities 04 and 
aj^^Z-Bgjy in order to establish scaled results for the maximum mass and central 

pressure and energy density. The value of /igo/^e// ^ quadratic solution of the 
equation for p = from Eq. (j44p . To compute numerical values for the mass, en- 
ergy density, pressure and chemical potential, it is then additionally necessary to fix 
Beff- However, this is not necessary to establish the quantity Smaxymax/^o which 
is independent of Eq — ^B^f / . 

We find that varying 04 for fixed nonzero (i'^'^ Ib]]^^^ has small effects on eigenval- 
ues hmax,Xmax, Umax and Smaxymax/^o- Likewise, there are similar small variations 
of these eigenvalues with varying /B^j^'^j: for fixed values of 04. In any event, the 
values of hmax,Xmax and Emaxylaaxl^o are largest for the case 04 = 1 and a-i = 0, 
as seen in Table [H Furthermore, the increase in the upper limit to the chemical 
potential resulting from 04 < 1 is approximately canceled by the ai contribution. 
Therefore, the upper limits to the central density and pressure established by the 
masssless quark bag model hold without regard to QCD corrections, pairing, or fi- 
nite strange quark mass. This result is confirmed plienomenologically, as illustrated 
in Fig. 2] which indicates that the limiting density curve for s — 1/3, to within a 
few percent, bounds pure quark stars and hybrid stars containing both quark and 
hadronic matter. As we proceed to show, the addition of a hadronic mantle and/or 
crust may permit densities slightly larger than that predicted by Eq. (|39]), just as 
they increased the densities for some hadronic neutron stars in some cases beyond 
the s — \ curve in Fig. |4l 



Simple Hybrid Stars 

The thermodynamic limits we have so far considered for quark matter apply to self- 
bound configurations with non-zero baryon density at their surface or those with 
a very thin hadronic crust. Hybrid stars consist of a quark core or a quark-hadron 
mixed phase together with a massive hadronic exterior. The presence of a mixed 
phase is possible only if the surface tension of quark matter is sufficiently large. We 
consider here the case where there is no mixed phase and the transition between 
quark and hadronic matter is a first-order transition at the point where /i and p in 
the two phases are equal, with values /it and pt , respectively. 
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We will employ a straightforward example. A reasonable approximation for 
hadronic pressures in the vicinity of nuclear matter densities for neutron star matter 
is p = Ke^, where iC ~ 1.1 • 10^** fm^ MeV~^. For h > ht we will use the general 
form for self-bound matter, Eq. (|18p . Requiring fi and p to match at ht, we find 

£ = (e''/2 - l)/K h<ht (45) 

e = eoA + e(''-''')*[£oA+(e'''/2-l)Vi^]As h > h. (46) 

In general, e and n will be discontinuous at ht, but we can require them to be 
continuous if we set 



1-- 



(47) 



This assumption does not have a significant effect on our results, and has the ad- 
vantage of reducing the number of free parameters in this model to a single one: 
ht, or, equivalently, pt- In the case that ht or pt tends to zero, we recover the self- 
bound cases, hmax increases slowly with pt, approximately as ln(l -\-pt/ps)- For the 
case Pt — lOps, integration of the TOV equation with this EOS gives values for the 
quantities in Table [2] in the rows denoted **. 

It is observed in both s = 1 and s = 1/3 cases with added shells and pt = lOps 
that hmax modestly increases by about 0.04, i.e., ^raax is increased by about 4%. 
Further increasing pt results in marginally larger values of hmax- Because the self- 
bound cores in both the s = 1 and s = 1/3 cases are stifFer than is realistic for 
hadronic and quark matter, respectively, which tends to unrealistically enhance the 
values of hmax, we conclude that the effects of adding a hadronic shell does not result 
in a significant increase in the central chemical potential, and the values established 
for the self-bound cases still hold approximately. However, values of £m,axVlnaxl 
decrease with the addition of a hadronic shell, so our model is not sufficiently general 
to explain the small violations of the s — \ and s ~ 1/3 curves seen in Fig. S) We 
now compare these results, in the case for quark stars, with those from published 
models. 



Published Models of Hybrid Stars 

^From many calculations in published literature it can be concluded that the den- 
sity at which exotica appear is uncertain and depends sensitively on the nature of 
nucleonic interactions up to (2 — 3) n^. In most cases, however, the EOS with exotica 
is softer than that without. Consequently, the maximum mass that a hybrid star 
can support is often smaller (with few exceptions) than that of a pure nucleonic 
star. In view of the observation of a 1.97 star, the natural question to ask is 
whether or not exotica are even permitted within such stars. 

At the current stage, it is difficult to rule out exotica as many model calculations 
indicate that a 2 Mq can indeed be supported by hybrid stars of different varieties. 
Indeed, Fig. |4] contains examples of hybrid quark stars exceeding this threshhold. 
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As another example, Fig. 5 of Ref. 3 shows kaon-condensed stars of 2 M© that 
would be compatible with the current observation. Ref. ,67 ^ which contrasts the 
MIT bag model with the Nambu Jona-Lasinio model for the treatment of quarks, 
presents results of hybrid stars containing nucleons and quarks, and nucleons, hy- 
perons and quarks (see Figs. 4 of Ref. [67)1 with maximum masses of ~ 2M0. The 
central densities of maximum mass configurations obtained there were in the range 
(0.7 — 0.75) fm~^ (see Fig. 5 of Ref. 'ST), which is compatible with the limiting curve 
s — 1/3 in Fig. |4l It wou ld be interesting to see whether or not Dyson- Schwinger 
-equation approaches to the quark matter EOS yield similar results. 

Recently, perturbative calculations of quark matter (without effects of pairing) 
including strong interaction effects to order and the effects of the strange quark 
mass have been performed by Kurkela et al. in Ref. |69l As these authors point out, 
their results become reliable at the 10% level only for quark chemical potentials 
exceeding 1 GeV, a value far larger than can exist in pure quark or quark hybrid stars 
as we have shown (the maximum is about 0.5 GeV). Their calculations nevertheless 
result in quark hybrid stars with masses approaching 2 Mq. 

These examples show that it is at present not possible to rule out any form 
of exotica with the 1.97 Mq observation. It must be stressed, however, that many 
of these approaches suffer from a great deal of model dependence (even in cases in 
which a substantially lower than a 2 Mq maximum mass star is predicted). However, 
the situation would change dramatically in the event of a confirmed observation of 
a 2.4 Mq star. Such an observation would seem to be incompatible with quark 
or quark hybrid stars including interactions, and would severely restrict models 
containing other forms of exotica. 

Gerry's Views on Hybrid Stars 

Prior to his fascination with kaon condensates, Gerry, together with Hans Bethe and 
Jerry Coopertseinl^, had maintained that perturbative QGD cannot be trusted at 
densities near rig where matter must be confined, and that any strangeness would 
appear in the form of Lambdas. To infer whether strangeness-rich quark matter can 
occur at high densities, they calculated the energy density of massless up, down 
and strange quarks to order as taking its density dependence as as = 2.2(n/nsY^^ . 
However, as shown in Ref. I66( their choice turned out to be independent of as 
once thermodynamic consistency was imposed. The result of their calculations was 
that the transition to quark matter occurred at densities larger than those found in 
neutron stars. 

A lesson lurks in Gerry's point of view. Most models of hyperon-filled stars have 
used hyperon-nucleon couplings based on simplistic quark counting in the non- 
perturbative regime. Analyses of laboratory experiments, however, indicate that at 
nuclear den sities t he A-nucleon potential is attractive, but the I]~-nucleon potential 
is repulsive l ^-^ l ^^ l. Knorren, Prakash and Ellis investigated the consequences of 
varying hyperon-nucleon couplings on the structural properties of neutron stars and 
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concluded that, if only Lambdas were to exist in neutron stars, their effects on the 
maximum mass would be minimal. If such is indeed the case, then two birds can 
be killed with one stone: first, a star with a 2-solar mass could be easily attained, 
and second, an avenue for a fast cooling could arise even if the direct Urea process 
with nucleons could not take place as even mi nisc ule amounts of Lambdas through 
their beta decays can cool the stars rapidly In a final example, based on a 
field-theoretical model for the hadrons and an NJL model for the quarks, Ref. [75] 
fits observed A and S~ hyper- nuclei and is able to construct hybrid stars of 2 Mq. 
Further investigations to fortify these ideas appear promising. 

5. Neutron Matter and Astrophysical Constraints 

It is also possible to place constraints on the maximum mass by combining the con- 
dition of causality with laboratory data concerning the symmetry energy, neutron 
skins on neutron-rich nuclei, and the properties of pure neutron matter as well as 
with astrophysical measurements of neutron star radii and masses. Recently, esti- 
mates of neutron matter properties have been obtained by first principle Green's 
Function Monte Carlo calculations and from chiral effective interactions 
These approaches yield similar results. Hebeler et al. showed that such calcula- 
tions restrict neutron star radii to be in the range of (11-12) km when the star's 
central density is at a density of l.lng, roughly the uppermost density for which 
the neutron matter calculations are reliable. With this constraint, it was possible to 
show that causality limits the maximum mass to about 2.9 Mq (see Fig. 3 of Ref. 

ESD. 

Steiner, Lattimer and Brown 1^ demonstrated that estimates of masses and radii 
of neutron stars in photospheric radius expansion X-ray bursters and also quiescent, 
cooling neutron stars limit the radii of 1.5 Mq stars to lie in the radius range (10.9- 
12.5) km. This radius range resulted from a Bayesian analysis of the observational 
data from three sources of both types. Encouragingly, t his rad ius range is also 
consistent with the results of neutron matter simulations l ^^ l ^^ l. With this mass- 
radius constraint, causality limits the maximum mass to about 2.4 M©, about 0.5 
Mq smaller than that inferred from the neutron matter results themselves. This 
upper limit too the maximum mass is due to the fact that the radius does not 
rapidly change with mass for a significant part of the mass-radius relation, as can 
be observed in Fig. [5J Thus, a restriction of the radius at a specific mass translates 
into an upper mass limit when the causality condition (i.e., M < i?c^/(2.9G)) is 
imposed. 

6. Conclusions 

Based on the observed largest well-measured mass, we have determined model- 
independent upper limits to the central energy density, central pressure and central 
chemical potential of neutron stars whether or not they are composed of hadrons. 
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quarks, or both. If matter undergoes a quark-hadron transition within the den- 
sity range permitted for neutron star interiors, these hmits on central density and 
chemical potential are substantially lowered. Indeed, if it can be shown that the 
maximum mass of neutron stars exceeds 2.5 Mq, the possibility of a quark-hadron 
transition in neutron stars is actually excluded, and the transition in cold matter 
would have to occur at an energy density in excess of SSs or about 1.2 GeV. 

We have also demonstrated the existence of an upper limit to the baryon chem- 
ical potential, independent of mass observations, of 2.1 GeV. If quark matter exists 
in the cores of neutron stars, the upper limit is likely not greater than 1.5 GeV. We 
also find the maximum binding energy or any neutron star is about 25% of the rest 
mass. 

The astronomical evidence currently favors the neutron star maximum mass 
to be between 2 Mq and 2.5 Mq. A Bayesian analysis of the observed neutron 
star mass distribution suggests that this distribution shows no indication of being 
influenced by an absolute upper limit imposed by general relativity and causality 
ISQl. In other words, this distribution shows no evidence of the existence of an abrupt 
upper mass cutoff; rather, the largest observed mass may simply be limited by the 
amount of mass it is reasonably possible for a neutron star to accrete in a binary. 
Supernova models generally imply that newly formed neutron stars are not more 
massive than about 1.5 Mq. In this case it is likely that the true maximum mass 
determined by the EOS is at least a few tenths of a solar mass larger than the 
largest observed mass. If so, the chemical potential, density and pressure limits we 
derived are somewhat conservative. In addition, our survey of EOSs with exotic 
matter indicates that stars containing significant proportions of this matter in the 
form of hypcrons, Bose condensates or quark matter would be almost impossible. 
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